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Abstract 

We consider maps between Riemannian manifolds in which the 
map is a stationary point of the nonhnear Hodge energy. The vari- 
ational equations of this functional form a quasilinear, nondiagonal, 
nonuniformly elliptic system which models certain kinds of compress- 
ible flow. Conditions are found under which singular sets of prescribed 
dimension cannot occur. Various degrees of smoothness are proven for 
the sonic limit, high-dimensional flow, and flow having nonzero vor- 
ticity. The gradient flow of solutions is estimated. Implications for 
other quasilinear field theories are suggested. 
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1 Introduction: nonlinear Hodge theory 



The original motivation for nonlinear Hodge theory was a conjecture of Bers 
on the existence of subsonic compressible flow having prescribed periods 
on a Riemannian manifold. In establishing Bers' conjecture for a station- 
ary, irrotational, polytropic flow/ L. M. and R. J. Sibner were led to more 
general questions about the properties of differential forms on Riemannian 
manifolds. In Ref. 2 they introduced a variational principle for the gen- 
eralized energy functional 

1 r rQ 

E^- / p{s)dsdM. (1) 



2 Jm Jo 

(See also Ref. 5, p. 221.) Here M is a smooth n-dimensional Riemannian 
manifold. Denote by a; a smooth section of the p*'* exterior power of the 
cotangent bundle on M; then Q{u!) is the square of the pointwise norm of cu 
on TM. The function p : R — > R+ is assumed to satisfy the conditions 

< ^ (Qp'iQ)) < oo (2) 



for Q < Qcrit, and 



If p = 1 and uj is the 1-form canonically associated to the velocity field of 
an adiabatic, isentropic compressible flow on M, then p is given explicitly by 
the formula 

/ ^ _ 1 \ V(7a-1) 

piQ) = (i - ^^Q) 

where 7„ > 1 is the adiabatic constant of the medium, and Qcrit — 2/(7a + l) 
is the square of the sonic flow velocity. 

The variational equations of E arc the nonlinear Hodge equations 

(5(p(QH = 0. (4) 

If the flow is irrotational then we have an additional equation 

d(jj = 0. (5) 

Here : ^ A^^^ is the (flat) exterior derivative on p-forms, with adjoint 
operator S : Ap ^ AP'^. 
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Applying the converse of the Poincarc Lemma to eq. (5), we find that, 
locally, there is a p — 1-form u such that du = lo. 

It has been observed^ that if {xi, are coordinates in R""*"^, u is 

a mapping of M into R""*"^ such that Xn^i = u{x^, a;„), p = 1, and 



then eq. (4) can be interpreted as the equation for a family of codimension- 
1 minimal hypersurfaces having gradients Vit. The critical value of Q is 
Q = oo. Of course E does not yield the area functional but rather an 
indefinite functional 



which differs from the area functional by an integration constant. 

\{ p = 2, n = 4, and uo denotes an electromagnetic field having elec- 
tromagnetic potential -u, then for p = 1 eqs. (4), (5) reduce to Maxwell's 
equations on M. If, however, we replace the standard model of electromag- 
netism with the Born-Infeld model, then we have (taking the energy to be 
positive-definite) 



where 6^ = mc^. The nonzero integration constant observed in (7) arises in 
(8) from independent physical arguments [c.f. eqs. (1.1) and (1.4) of Ref. 6]. 
Normalizing so that 6^ = 1/2, the energy functional (8) becomes identical to 
the functional (7), and we can choose p as in eq. (6) in order to write the 
variational equations of this common energy functional in the form of eqs. 
(4) . Thus the Born-Infeld model fits naturally into nonlinear Hodge theory as 
an application for 2-forms. (The equations for a nonparametric codimension- 
1 minimal surface also have a place in the original gas dynamics context of 
nonlinear Hodge theory, as the Chaplygin approximation of a compressible 
fiow.^) 

Because the bundle T*M is fiat, any connection defined on it will have 
trivial Lie bracket. For this reason, in comparison to the examples that 
follow, we call the foregoing cases abelian. In particular, in the example 
of electromagnetism the vector potential u is identified with a connection 
1-form on a bundle over M having abelian structure group U{1). 



p(Q) = (i + Q)-'/', 



(6) 




(7) 
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Suppose we replace the energy functionals (1) and (7) by the functional 



where 



p{s) = {h + s) 



—a 



Here h and a are nonnegative parameters. Let X be a vector bundle over 
M having compact structure group G. Define Q = {F, F) , where ( , ) is an 
inner product on the fibers of the bundle adX <S) hP{T*M),p > 1. 

Case 1: Let p = 2 and let A be a connection 1-form on the fibers of 
X; choose F to be the curvature 2-form Fa corresponding to A. If o; = 0, 
G = SO{n), and n = 4, then Eq is the Yang-Mills functional. 

Case 2: Ifp = 2,Q! = l/2,Xis the bundle of orthonormal frames on 
M, G is the Lorentz group 3), and n — A, then Eq is formally analogous 
to a (torsion-free) gravitational action functional for the curvature 2-form F. 

Case 3: If p = 1, a = 1/2, z = A{x, y), where A is the graph of a surface 
E in R^, F = gradA, X = T*M, and n = 2, then Ei is the energy functional 
for a nonparametric minimal surface, as discussed above. 

Case 4' If in the last example is replaced by the Minkowski space 
R'*'-^ and if p = n — 1, then Ei is closely related to an energy functional for 
maximal space- like hypersur faces. ^ 

With the exception of the last two, these examples are characterized by 
a nonvanishing Lie bracket in A due to a nonabelian structure group for X. 
Thus in general these cases are nonabelian. If Da represents the exterior 
covariant derivative associated to a connection 1-form A on ad{X) and if D\ 
is the formal adjoint of Da, then the variational equations for Eh can be 
written®" 



The first equation represents a nonabelian version of eq. (4) for curvature 
2-forms, and the second equation replaces eq. (5) by the second Bianchi 
identity. 

An intermediate place between the abelian and nonabelian nonlinear 
Hodge theories is occupied by nonlinear Hodge maps. These are maps u 
between Riemannian manifolds such that w is a critical point of the non- 
linear Hodge energy (1). In this case the geometry of the target space is 




DaFa = 



0. 



(9) 
(10) 
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enriched in comparison to the abehan case but does not have the nontrivial 
Lie group structure of the nonabchan case; the target space is independent of 
the base space but is not a curved bundle. In the context of fluid dynamics 
or the rotation of a nonrigid body, these maps represent flows on a Rieman- 
nian manifold M (typically, a domain of R") for which the flow potential is 
constrained to lie on a possibly different Riemannian manifold N. 

Nonlinear Hodge theory can be viewed as an attempt to extend to the 
quasilinear field equations of classical physics the unified geometric treatment 
given linear field equations by the theory of Hodge and Kodaira. The case 
p = 1 for the abelian equations (4), (5) reduces to the continuity equation for 
an incompressible flow {p — 1) or the field equations for electromagnetism 
(p = 2). For the nonabelian equations (9), (10) the case p = 1 reduces to 
the Yang-Mills equations. In the intermediate case considered in the sequel 
the case p = I reduces to the equations for harmonic maps (nonlinear sigma- 
models). In distinction to the approach in Refs. 1-4, our concern is with the 
geometry of the target space rather than the geometry of the domain. 

In the following we denote by C generic positive constants, which may 
depend on dimension and which may change in value from line to line. We 
employ the summation convention for repeated indices. 



2 The variational equations 

Consider a map u : M ^ N taking a Riemannian manifold (M, 7) into a 
Riemannian manifold (N.g). We are interested in maps which are critical 
points of the nonlinear Hodge energy (1). Here and throughout we denote 
by a; = (x^, x") a coordinate chart on the manifold M having metric 
tensor 7a^(a;), and we denote by u — {u^, a coordinate chart on the 

manifold N having metric tensor gij{u). (We assume for the moment that 
the image of u lies in a coordinate chart.) The nonlinear Hodge energy 
assumes the form 

I r fQ 
E{u) = - / / p{s)ds^dx, 

Z JM Jo 

where 

Q = rn.)9« («(.)) 

We have, by the Leibniz rule, 

^E{u + t^P)\t=o = ]- I p{Q)^Q{u + t^P)\t=o^dx, (11) 
dt z JM dt 
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for arbitrary ip G C^{M). The construction of the test functions ip is not en- 
tirely straightforward. Use must be made of the Nash Embedding Theorem 
to embed in a higher-dimensional euclidean space. One then employs a 
nearest point projection 11 of a suitable euclidean neighborhood 0{N) onto 
N. If t is small enough and N isa,C^ submanifold, the variations IIo (u + ttjj) 
will be constrained to lie on A^, where now ip : M ^ O. A discussion is 
given in Section 1 of Ref. 11 for the special case of harmonic maps. 

Carrying out the indicated operation on the right-hand side of (11), we 
obtain 



If , , . d , ..... ,i.du^ du^ 



But 



Jm 



I 



,^p(g)7"^(x)^ (^.- {u{x)))^P^^^V^dx. (13) 

Substituting eq. (13) into eq. (12) and taking into account that t/j has 
compact support in M we obtain, for 77* = Qij'ip'' , the formula 



2 Jm^^^'^ ^ ^ ' \dx^ ^ dx^ dxo ) dx- dxP^ ^^"^^^ 

Applying the definition of affine connection, we conclude that stationary 
maps satisfy the system 

' ^ ^(Q);^7"^^(a:)|^| + P(g)7"^(a:)r;,(.)|^^ = 0, (14) 
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for i = 1, m. 

Wc can also approach the variational equations for nonlinear Hodge maps 
from an intrinsic point of view, defining the nonlinear Hodge tension field r 
by the formula 

T = trace Wcov {p{Q)^^) , 

where Vcov denotes the covariant derivative in the bundle T*M ® u^^TN. 
li N — H'^, then the equation r — reduces to the conventional nonlinear 
Hodge equation (4) for uj = du [which implies eq. (5)]. In particular, 
if H o (-u + tip) are the variations described earlier, then we can write the 
equations for a weak stationary point in the form 

I I p{Q) [{Vu, V^) + V^u {DU{u) [Vu] , DU{u) [Vm])^] dM = 0, (15) 
where 

I {Dn{u) [Wu] , Dn{u) [Vu])^| < C IVI \u\ I Vk|' ; (16) 

Wu = grad u; V.,pu { , is the function on M whose value at x is the 
covariant derivative for of the metric ( , )^ in the direction iIju{x). (See 
the discussion leading to inequality (1.2) in Ref. 12; that paper considers the 
regularity of energy minimizing maps for p{Q) = Q^-) Here we use the fact 
that for It e A*^, I a; I is the norm of the gradient of u as well as the norm of 
its differential. 

The harmonic map density satisfies 

( \ =- ocP I ^ W \ 

^[U ) harmonic — ^1 ^ ' dx^ / ~ 2 ^ ' ' 1^*-^®""^^-^ ' 

That is, in the harmonic map case the energy density is the trace of the 
pullback, via the map u, of the metric tensor g[u) on A^. In the case 
of nonlinear Hodge maps the situation is a little more complicated, as the 
energy density is the integral 



F(a;)= 1'^ p{s)ds, 
Jo 



which may not be a quadratic form if p 7^ 1. Moreover, the nonlinear Hodge 
density need not scale like a metric tensor. Thus the geometry of harmonic 
maps is more transparent than the geometry of nonlinear Hodge maps. 
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Proposition 1 In order for weak solutions uj = du of the equations r — to 
exist locally on M it is sufficient that N be there exist a positive constant 
K < oo for which 

<p{Q)+2Qp'{Q)<K. (17) 

Proof: The argument follows Sec. 1 of Ref. 13. Define F((j) as in the 
preceding paragraph. Then (17) implies that 

>0 



and there exist finite positive constants ko and ki such that 

koQ < — - — (jJaOJp < kiQ. 

OUJf3OL0a 

Moreover, there exist finite positive constants k2 and ks such that 

k2Q < F {Q{uj)) < hQ. 

Thus the energy functional E is convex, bounded above and below, and 
lower semicontinuous with respect to weak LF' convergence. This completes 
the proof of the proposition. 

We will use in several contexts the following pointwise inequality for 
smooth solutions. 

Theorem 2 Let u : M ^ N be a stationary point of the nonlinear Hodge 
energy on M, where M is a compact, n- dimensional C°° Riemannian mani- 
fold, n > 2, and N is a compact m- dimensional C°° Riemannian manifold. 
Then the scalar Q = IVu] satisfies an inequality of the form 

L{Q) + C{Q + 1)Q>0, 

where the second-order operator L is elliptic whenever condition (17) is sat- 
isfied, and the constant C depends on the Ricci curvature of M and the 
Riemann curvature of N. 

Proof. Denote by a subscripted x"^ differentiation in the direction of the 
^th coordinate. Differentiation of the metric tensor and Christoffel symbols 
in the direction of an index is indicated by a comma preceding the subscripted 
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index. Choose geodesic normal coordinates at the points x & M and u{x) e 
N. At these points 

r^x) = S^^ix); g,,{u) = 5,,{u)- T%{x) = T1^{u) = 0. 

As in the preceding, greek indices are used for coordinates on M, and latin 
indices, for coordinates on A^. 
Write eq. (14) in the form 



{p (Q) + p' (Q) Q,.ul, - p (Q) [uir,Tl^{x) - rl{u)ul.uU 
Differentiating (18) with respect to x^ and letting a = (3 = a yields 

P (Q) '^X^X'^X^ — ~P' (Q) Qx^ul-a^cT — p' (Q) Qx'^'^x 

+ P (Q) [uivT'l^^.ix) - rl^^p{u)ulsu^^.u^ 
Now compute (c./. Sec. 3.2 of Ref. 14) 

Ae{u)=[r''9ij{^)p{Q)uUui 
^''^gij/kiu)ui.u^^aP (Q) 4«</3 + [t^^^^jP' (Q) Qx'^uUu 

I'^^gijp iQ) K'^x'^K^]^. + ^"^9ijP (Q) <c.x'^Kex''+ 

7 

Here 



}^0. 
(18) 



(19) 



4/3 



+ 



T4 = 



^p iQ) Qx- 



(20) 



(21) 



Subtract this term from the left-hand side of eq. (20). Applying eq. (19) 
with e — a, T] — P, we have 



where 



Ti + T2 + Te + Tg = 

HQ) = {p' iQ) Ka [Qx-ui^ - Qx'^ui^] , 
9 



(22) 



R^p is the Ricci curvature of M, and -R^^-^ is the Riemann curvature of N. 
The last term in (22) results from applying the product rule to the quantity 
—u^x-^ [p' (Q) Qx'^u'-^Ax" ■ write 



and 



|T*M(8iu-irW 



If p' {Q) < 0, then (17) and the Schwarz inequality imply that 



If p' (Q) > 0, then 



75 + TV > 7^5 = p(<5) I Va;|' > C | Va;| 



(23) 



(24) 



[See the argument leading to (26), below.] In either case we obtain from 
expressions (20)-(24), defining 

L{Q) = Aeiu) - T4 - A{Q) 



p{Q) + Qp'{Q) 



qA -A(Q), 



the inequality 

HQ) > p{Q) {R^puUui, - R^^,v},.uUul.ui.) + C . 
We can write 

L{Q) + C$g > 0, (25) 

where $ = Q + 1- 

We now show that the operator L is eUiptic whenever condition (17) is 
satisfied. If p'{Q) < 0, then 

— p' [Q) Ul^a Qx'^'^x" — Qx"'^x<^ Qx" — 

\P (Q)!"^!;" ^Qx'^ui.a — Qx^U^x" Qx" — 
\P' {Q)\ {Q |VQ|' - > 

10 



p'(Q)l(Q|vg| 
1 



|p'(Q)l(Q|vg|^-- 



Qx'^Qx" ) ^ 



0, 



where in the final step we have applied Young's inequality. Thus in this 

case 

":P(<9) + Qp' (Q)] |vg|' - p' (Q) uU IQx^uI^ - Qx^ > 



\VQf>-K-'\VQ\' 



If p'iQ) > 0, then p{Q) > p{Q^^) = p(0). At Q = 0, 

\p{Q) + Qp' {Q) = \p- 

Condition (17) then implies that p is bounded below away from zero. Again 
using Young's inequality, 



(26) 



i^p{Q) + Qp' {Q) 



> 



p{Q) + Qp' {Q) 



|vg|^ + p' (Q) v:^cAj:i,.Qx^Qx- > 



-p(Q) + Qp' (Q)J I VQI' - Qp' (Q) \VQ\' = -p \VQ\' > C \VQ\' , 

as p is bounded below away from zero. 

Thus for either sign of p'{Q), condition (17) implies that there is a con- 
stant mi > such that 



^p{Q) + Qp' (Q) 



and we can write 



\WQf - p' (Q) [Qx^uU - Qx^ui.] Q^. > mi | Vg|' 



HQ) 



d 



dQ' 



for a matrix a'-' satisfying 



mi lel' < a'^Uj < m, \^\' 



(27) 



(See p. 106 of Ref. 3 or Proposition 1.1 of Ref. 4 for a different proof of 
ellipticity.) This completes the proof of Theorem 2. 

The contribution of the target geometry to the nonlinearity of inequality 
(25) significantly exceeds that of the geometry of the base manifold, which 
in the sequel we generally take to be euclidean. 
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3 Apparent singular sets of prescribed dimen- 
sion 

In general we expect that finite-energy weak solutions to cqs. (15) may be 
singular, as singularities occur even in the case p = 1. It is natural to seek 
extra conditions under which solutions are actually smooth. 

Theorem 3 Let u : Q ^ N be a C"^ stationary point of the nonlinear Hodge 
energy on Vt/H, where Vt is a domain of R", n > 2; N is a compact Tri- 
dimensional C°° Riemannian manifold, m < n; T, GG B GG fl is a compact 
singular set, completely contained in a sufficiently small n-disc B, which is 
itself completely contained in Let conditions (16) and (17) hold. Ifn > A, 
let 2n/ {n — 2)<fi<n, where p is the codimension ofH, and let d,u G L^{B). 
// n = 3, 4, let du e n where /S = {p - e) / {p ~ 2 - e) for 

2 < p < n, e > 0, and ^ < Qq < q. Then du is Holder continuous in fl. 

Remarks. That an condition is necessary even for p = 1 can be seen 
by considering the equator map.^^ In Theorem 3 we assume neither that 
the map u minimizes energy nor that the energy is small. In distinction to 
the harmonic map case, we do not obtain higher regularity from the Holder 
continuity of du, as the system (15) is not diagonal in it principal part. In 
the following proof wc assume that p — n; the extension of the proof to lower 
values of p can be effected by arguments given in Sec. 3 of Ref. 27. 

Proof. Initially let n exceed 4. Integrate inequality (25) against a 
nonnegative test function ( e C^{B) given by 

C = (#)'s(g), 

where B — Br^xq) is an n-disc, of radius < i?, centered at a point Xq e 
O; assume that B completely contains S and is completely contained in 

^ ^'li^ix) = OVx in a neighborhood of S;?] G C^{B') where B' GG B; 
S(Q) = H{Q)H'{Q), where H{Q) = Hk{Q) is the following variant of Serrin's 
test function^^: 



g[n/(n-2)]«n/4 forQ<Q<l, 

. g(M-2-.)/2y"/("-2)r"/2(/^--) _ ^^[n/(n-2)]'^n/4 



IJ.-2-e 



for Q>e. 
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Iterate the following sequence of elliptic estimates, taking successively u e 
for a{K) = n[n/{n - 2)]^, k = 0, 1, ... . For all k < oo, 



a''{u)diQ-2{7]ilj)dj{7]ilj)E{Q)*l 

+ [ a'^{u) {r]^f E'{Q)diQdjQ * 1 < / $g E{Q) * 1. 
Jb' jb' 

This inequality can be rewritten in the short-hand form 

h + h<h: (28) 

the integrals of which we estimate individually. Because exceeds 2n/ {n — 2) 
we have 

E'{Q)>C{H'{Q))\ (29) 

Also 

Inequality (29) implies, by ellipticity, 

h = a'^{u) {rii^f E'{Q)diQdjQ * 1 > 

C(mi) / (r7V)'(i/'(Q))'|Vg|%l = 

JB' 

C I {r)'ipf\WHf*l = i2i. (31) 
Jb' 

Young's inequality implies 

h = J^^ a''{u)diQ ■ 2 {fjiP) (rjiP) H{Q)H\Q) * 1 = 

2 (a'^(M) (#) (a,i/)) a,- (#) * 1 > 

-m2 (#)^ |Vif|^ * 1 + C{e) j^^ |V (#)|^ //^ * 

= - (ill + ii2) . (32) 
Using inequality (30) and the Sobolev inequality, we obtain 

h= [ ^Q(ri^fE{Q)*l<(^yj[ $(#)^i/2*l< 
JB' \n — 2/ 4 JB' 

13 



< c \mL/2 { J^, [I V + m'] H'*i + M)' I vi/|' * 1} 

= ^31+^32- (33) 

For sufficiently small B'we have 

0<i2i-(i32 + iii) <C'(ii2 + i3i). (34) 

There exists ( c.f. Ref. 17, Lemma 2 and p. 73) a sequence of functions 
such that: 

a) C e [0, 1] Vz/; 

6^ = 1 in a neighborhood of S, Vi/; 

— > a.e. as — > 00; 
c?; V^^ ^ in L^"-^ as 1/ ^ 00. 

Apply the product rule to the squared H^'^ norm in i^i letting ^p = ?,/-'^ = 
1— Observing that the cross terms in (Vf]) ip and (Vip) f] can be absorbed 
into the other terms by applying Young's inequality, we estimate 



lim / rj^lW^ljf *1 < lim C(£) [ IVz/'J^Q'^^^)^ 



2 * 1 



< Jim C(£) ||V^J|^_, \\u\\i$,^'-'-'^'^'-'^ = 0. (35) 

Having shown that the integral on the left in (35) is zero for every value 
of we can now let £ tend to infinity. We obtain via Fatou's Lemma the 
inequality 

1 V 1% 1 < j^^ I vr;|=^ g"W/2 * 1. 

Thus Q'^^'^^/^ is in H^'^ on some smaller disc on which 77 = 1. But then, 
because u is assumed to be away from the singularity and S is com- 
pact, must be in H^'^ on the larger disc as well. Apply the Sobolev 
inequality to conclude that u is now in the space L°''^'^^^\B) . Because the 
sequence {n/{n — 2)}'^ obviously diverges, we conclude after a finite number 
of iterations of this argument that Q'^ is in H^'^{B) for any positive value of 
c. A final application of the Sobolev inequality implies that uu £ L^{B) for 
all s < 00 and for any small B <Z<Z fl, provided n exceeds 4. 
Now let n = 3 or 4. Dcfine^^ 
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Q1' forO<Q< 



1 r„/ 



90 



[q'ii'-QoQqo _ (^q^ _ q'^ for Q > £, 



where q' = [n/ {n — 2)]'^ q. Arguing as in the higher- dimensional case we 
obtain, using the Sobolev inequahty, Q E L'^i ^/^^'■^'^^(^B). Repeating the 
argument for = 0, 1, . . . , we obtain that uu € L^{B) for all s < oo when n 
is 3 or 4. 

Now let n be an arbitrary integer greater than 2. Again let ip^, — 1 — 
where satisfies properties a)-d) above, and let rj e C^{B') as before. If 
C = V'z., then 

/ {dC, p{Q)uj) *1= f (C, piQ)biuj)) * 1, (36) 
with b given by (15). We have 

/ {dC, p{Q)u) * 1 = / U (d^,) , p{Q)u;) * 1 

JBu J Br \ / 

+ j^^ {i^,d{ri'),p{Q)u;)*l, 
where as u tends to infinity, ip^, tends to 1 a.e. and 

/ {v' {di^;} , p{Q)u:) * 1 < C{K) - 0. 

J Bji 

Choosing rj^^x) to equal 1 for x G we find from (36) that uj is locally a 

weak solution in all of fi. 

Let the map tp : Br{xo) — > R*^ satisfy for sufficiently large k the 
boundary-value problem 

S {p{Q{d(p))dip) = in Br{xq), 
(fi^ — on dB, 

where the subscripted ■»? denotes the tangential component of the map in 
coordinates (r, 'di, ...i^^-i). The existence of a C^'" solution ip to this problem 
is well known. ^ Moreover, if {dip)ji^xo denotes the mean value of the 1-form 
dip on Br{xq), then dip satisfies a Campanato estimate^^ 

[ \dip-{dip)R,,f*l<CR'+^^» 
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for some number 7^ e (0, 1]. Then u — (p is an admissible test function, and 

/ {d{u- (p) , [p{Q{du))du - p{Q{d(p))d(p]) * 1 

= / {{u-ip),p{Q)b{u,Du))*l (37) 

with b given by (15). Apply to identity (37) Sibner's mean- value formula 
(Lemma 1.1 of Ref. 13), which asserts for the unconstrained case that 



G"(e, /) - G"(r/, h) = A'^^ iff, - hp) + - 



(38) 



where 



dF 



A"'!^ is a positive-definite matrix, and 

H^\<C{\f{x)\ + \hix)\). 

Here F is the function used in the proof of Proposition 1. Equation (38) 
extends immediately to our case, as we can estimate the derivative of the 
metric g on N hy 



dg 



dx 



< 



dg 



du 



du 



dx 



<C\u\ 



[c.f. inequality (1.3c) of Ref. 13; the g in the above inequalities is not the 
same object as the g in Ref. 13, which corresponds to our 7]. 

In formula (38) choose ^ — x, r] — f — u, and h — d(p. We obtain, 
using (16), 

/ \d{u — ip)f *1 < 

Cil {\u;\ + \d<p\)\x\*l + [ \u-<p\p{Q)\u\Q*l] . (39) 
We can find a number s sufficiently large so that 

/ {\^JJ\ + \d^p\)\x\*\< 

JBr{xo) 



R 



\X 



s/(s— 1) I in— 1 



(s-iys 



I n — 1 7 I I 

x\ d \x\ 
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< C{s, = CEP, (40) 

where rj > n whenever s > n. Also, Young's inequahty yields 

/ I"" — p{Q) \u\Q * I < 

JBr{xo) 

R"" [ lu-ipfluf^l + R^f QV(g)*i< 

JBr(xo) JBr{xq) 

R-'' f \u-ip\'\uf*i + c iWQi , WpU r-^'-'^''^^ 

JBii{xo) 

for a constant u to be chosen and s so large that us > n. We have 

R-" [ \u-ipf\uf*l< 

JBr{xo) 

/ \ (n-2)/n / ^ ^ 2/n 

R 



\JBr{xo) J \JBr{xo) , 

<R-''Csf \V{u-ip)f*li f l^r*! , (41) 

JBr{xo) \JBr{xo) ) 



where Cs is Sobolev's constant. The hypothesis on d\i now implies by the 
Sobolev Theorem (and a trivial application of the Gaffney-Garding inequal- 
ity) that for any £ > we have 



Br{xo) 



1 


7 M 




UBiiixo) 



\n/{n+e) / \e/{n+e) 

n+e r^-^dr / r"-Mr < CR^ 

J \JBr{xo) J 

for X — ne/ (n+e) . Because of the high W space in which u sits we have some 
flexibihty: choosing either £, s, or u so that u < X allows us to subtract the 
right-hand side of inequality (41) from the left-hand side of inequality (39). 
Because the mean value minimizes variance over all location parameters, we 
find that 

/ ^ \uj - {uj)R,^,f *1 < \uj - {dip)R,^,f *1 

JBr{xq) JBr{xq) 

< \uj-d(pf*l+ \d(f - {d(f)R^^Qf *1 

JBr{xo) JBr{xq) 

Choosing xq so that E CC Sij(xo), and observing that if = n, this argument 
works for any smaller positive value of R, completes the proof. 
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4 The sonic limit 



Denote by 7^ a closed 1-form having prescribed periods. We add to eqs. 
(4), (5) the homology condition^ that a; — 7^ be an exact form. Denote by 
M a smooth, compact Riemannian manifold and consider a family of maps 
Ut : M ^ N into a smooth, compact Riemannian manifold A^. We further 
assume that for each t : < t < tcrit, = dut is a weak minimizer of 
the nonlinear Hodge energy on M in the following sense: condition (17) 
is satisfied, eqs. (14) are weakly satisfied by the vector field canonically 
associated to Ut, tOt — t^i is an exact form in LF'{M), and for all other 1- 
forms a e L^{M) such that a — t^i is exact, the inequality 

/ / p{s)dsdM< / p{s)dsdM 
Jm Jo Jm Jo 

is satisfied. Borrowing the terminology of fiuid dynamics^ we call weak 
solutions cut, t e [0, tcrit), subsonic. The question is whether such solutions 
converge, as t tends to tcrit, to sonic solutions having velocity Qcrit- EUiptic- 
ity degenerates in the limit as Q tends to Qcrit [cf- eq. (3) of Section 1]. In 
this limit condition (17) fails and is replaced by conditions (2), (3). In the 
following theorem we replace M by a euclidean domain; but see the remarks 
at the end of this section. 

Theorem 4 Assume the hypotheses of the preceding paragraph. That is, let 
Ut : fl ^ N denote a family of maps between a smooth, compact domain Q of 
R" and a coordinate chart on a smooth, compact m-dimensional Riemannian 
manifold N, m < n,where < t < tcrit- Assume that the 1-forms cut — duf 
weakly minimize the nonlinear Hodge energy on Q over a cohomology class. 
In particular, let the homology condition of the above paragraph be satisfied 
for a fixed 1-form 7^. Assume that the function p satisfies (2), (3) and 
that 

Q<C [ [ p{s) ds *1 Vg < Qcrit (42) 

Jq Jo 

for constant c. Then as t tends to tcrit, 

lim max Q {cutix)) ^ Qcrit- (43) 

t^tcrit int fl 

The conclusion of Theorem 4 implies that Ut depends continuously on 
t in the topology of uniform convergence. This eventually implies Holder 
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continuity for weak minimizers at the elliptic degeneracy represented by (3); 
see Corollary 5. 



Proof. The proof is similar to that of Theorem 4.8 of Ref. 13. Denote by 
{ti,} a nonnegative sequence of points in [0, tcru) converging to a limit point. 

We want to establish a sequence of inequalities satisfied by any subsonic 
minimizer ujt^ = uju- Because lUi, minimizes energy over a cohomology class 
we have 

c / p{s)ds * 1 



in Jo 

< c / / p{s)ds * 1 < C \\K\\L2(^n) , 

where is a harmonic form such that h„—t„ji is exact. This gives a uniform 
bound in L°° on the sequence {cUu}- Now we proceed as in the concluding 
arguments in the proof of Theorem 3, comparing cui, to a solution (p of the 
euclidean nonlinear Hodge equations on Q. Conditions (2) and (42) imply 
that u — (f is an admissible test function. The continuity estimates, starting 
with formula (37), are also uniform, as the highest bound imposed on by 
these inequalities is in L""*"^. For example, we can replace inequality (40) by 
the estimate 

/ {\uj\ + \d(p\) \x\ * 1 < 

/ R \ {n-l)/n 

c(iieiL/2+ wu) (/„ w"/("""|xr'<i|x|j 

<c[r^ +\\d^\\^Aj^ \x\'"U\x\\ ^^K'KCW 

for some rj > n. Thus the hypotheses of Theorem 4 imply the concluding 
Holder estimate of Theorem 3, from which we obtain equicontinuity for the 
sequence {o;^}. Now the Arzela-Ascoli Theorem guarantees uniform con- 
vergence of a subsequence to a 1-form satisfying both the equations and the 
homology condition. This completes the proof of Theorem 4. 

Corollary 5 Let the hypotheses of Theorem 4 be satisfied. Then cu is Holder 
continuous in the interior ofD,. 
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Proof. Theorem 4 is the crucial ingredient in the technique of Shiffman 
regularization, described in the Appendix to Ref. 13. This technique is 
sufficient to estabhsh the Holder continuity of uj and prove the corollary. 

In Ref. 10 a comparison argument similar to (37)-(41) was constructed for 
solutions of eqs. (9), (10). There we chose an exponential gauge at the origin 
of coordinates in a euclidean ball -B_r(0) in order to compare solutions of (9), 
(10) with euclidean solutions of (4), (5). It was shown that the difference of 
the two solutions is small in a high Campanato space. It was then necessary 
to show that the gauge transformation to an exponential gauge preserves the 
Campanato estimate; this allowed us to extend the comparison outside of 
-Bij(O) and apply a covering argument. The argument of Ref. 10 provides 
a guide for extending the results of this section to maps of a Riemannian 
manifold M. The analogy of an exponential gauge is a choice of geodesic 
normal coordinates in a local coordinate chart. The arguments of Ref. 13 
imply that the difference between a comparison map (p, taking a euclidean 
ball into R'^, and a comparison map ip', taking a Riemannian ball into R^, is 
itself small in a high Campanato space. This is the analogy of our estimates 
of the gauge transformations in Ref. 10. Now we can extend the local 
estimate to all of M by a covering argument. Although in principle this 
method could be used to extend Theorem 3 to a Riemannian domain as well, 
in that case no covering argument is needed because S is assumed to be 
small. 



5 An application to harmonic maps 

We now consider the special case in which E is a point, p is constant, and n 
exceeds 4. The following result is a special case of a theorem which Liao^° 
proved by quite different methods. 

Theorem 6 (Liao). Let u : Q ^ N be a C'^ stationary point of the non- 
linear Hodge energy with p = I on fl — {po}, where fl is a domain of R", 
n > A; N is a compact m- dimensional C°° Riemannian manifold, m < n; 
Pq & is a point. If Q — \du\'^ satisfies the growth condition 

7o 



Q{x) < 



|2 



F - Pol 

for X e Bji(po), where Br(j)q) is an n-disc of radius R centered at po and Jq 
is a sufficiently small positive constant, then du is Holder continuous on Q. 
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Proof. The growth condition guarantees du G L^{B) VP < n. The idea 
of the proof is to show that du G L'^B) and apply Theorem 3. Without 
loss of generality we take po to lie at the origin of coordinates in R"^. 

Let^i ^{x) = C{x)'iIj{x), where x G Br{0) - {0}, 



ip{x) — \x 



i—n 



and C is chosen so that C,{x) = 1 \l2e < \x\ < i?/2, and C,{x) = if |a;| < £ 
or \x\ > R. We can find C, satisfying the additional conditions that 

ivci < f 

and 

Because L is a divergence-form operator and has compact support in Br, 
inequality (25) implies that 

-/ (A,OQ*i = -/ e(Ag)*i 

<C f * 1, (44) 

where is the Laplacian in radial coordinates. We have 

where 

AV' = 2(4-n)|x|'"". 
We can write inequality (44) in the form 

/ i^CQ(^-CQ\*l<2 f IVCIIVV'IQ*! 
Jbr{o) \ ip / Jbr{o) 



+ |AC|VQ*1- (45) 

Jbr{o) 

We are interested in the behavior of this inequality as the constant e in the 
trapezoidal function tends to zero. Write (45) in the form 

ii <2i2 + is. 
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Because in Bji{0) — {0} we have Q < 7o ^ , integration in radial coordi- 
nates yields 

C f 

i2 < J J^d\x\+C{R), 

where 

r = {x\e< \x\ < 2e}. 
Integral 12 is obviously finite as £ tends to zero. Similarly, 

C r 

«3 < ^ J^\x\d\x\+C{R), 
which is also finite for every e. Finally, 

The quantity inside the largest parentheses on the right is positive provided 
7q is sufficiently small. In this case 

limii > C Q * 1. 

But also, 

/ g-/2 * 1 = / Q (g(-2)/2) , 1 



7 Q 



Taken together, these inequalities imply that uj lies in the space in a 
neighborhood of the singularity. The hypotheses of Theorem 3 being satis- 
fied, we conclude that a; is Holder continuous, which completes the proof of 
Theorem 6. 



6 Rotational fields 

In this section we study systems of the form 

5(p(gV) = o, (46) 

duj /\uj, (47) 
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where to e (T*M) for p > 1; f G t^{T*M)\ M is an n-dimensional 
Riemannian manifold; Q = {uj,uj) = *(uj A *uj); * : ^ jg tJ^e Hodge 

involution; p : R — > R+ is a function satisfying the condition^ 

K-\Q + ky <p(Q) + 2Qp'{Q) <K(Q + ky (48) 

for some positive constant K and nonnegative constants k, q. 

If V = (or if p = 1 and v = lu), then condition (47) degenerates to 
condition (5). If a; e A^(T*M) is the 1-form canonically associated to the 
velocity field of an n-dimensional fluid, then condition (5) guarantees that 
the flow is irrotational: no circulation exists about any curve homologous to 
zero. 

If G A^(T*M), then condition (47) only guarantees, via the Frobenius 
Theorem, that cu — idu locally; a potential exists only along the hyper- 
surfaces i = constant, and circulation about topologically trivial points is 
excluded only along these hypersurfaccs. (For the extension of this result to 
exterior products of 1-forms, see, e.g., Ref. 22, Sec. 4-3.) Equations (4), (5) 
can be used to prescribe a cohomology class for solutions as in Sec. 4, but 
eqs. (46) , (47) will only prescribe a closed ideal. 

We have as an immediate consequence of (47) the condition 

du;Au} = 0. (49) 

If (J denotes tangential velocity of a rigid rotor (p = p{x) only), cq. (49) 
corresponds in three euclidean dimensions to the fact that the direction of 
V X a; is perpendicular to the plane of rotation. Condition (49) also arises 
in thermodynamics. ^^'^^ 

As in preceding sections, we replace M by a euclidean domain in proving 
the technical results. In the general case, the curvature of M enters in a 
predictable way.^^ 

Theorem 7 Let u!,v smoothly satisfy eqs. (46), (47) on a bounded, open 
domain Q C R". Assume condition (48). Then the scalar Q — *{u! A *a;) 
satisfies the elliptic inequality 

LUQ) + C{Q + ky{\Vv\ + \v\^)Q > 0, (50) 

where is a divergence-form operator which is uniformly elliptic for k > 0. 
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Proof. We have (Ref. 5, (1.5)-(1.7)) 

{00, A {p{Q)Lo)) = d; {lo, d, (p(Q)cu)) - {d^Lu, d; {p{Q)uj)) 
= ^H{Q) - [p{Q) {diu, diu) + p'iQ) {diu, u) diQ] , (51) 

where 

AH{Q)^d, 



^p{Q) + Qp'{Q))diQ 
e Q. Observe that H is defined so that 

H'{Q)^^p{Q) + Qp'{Q). 

Just as in the derivation of inequahty (25), we have 

p'iQ) {diu, u) diQ = 2p'{Q) {diu, uf . (52) 

i 

If p'{Q) > 0, then (52) imphes that 

p{Q) {diLo, diUj) + p'{Q) {diUj, u) d^Q > 

p{Q)\Vijf>K-\Q + ky\Vu\\ (53) 
In (53) we have used the inequahty 

p{Q)>K-\Q + ky, (54) 

which follows from (48) (with a possibly larger constant K). If p'{Q) < 0, 
then (52) and the Schwarz inequality imply, just as in the derivation of 
inequality (25), the inequality 

p{Q) {diUj, diu) + p'{Q) {diUj, u) diQ > 

p{Q) \Vuj\^ + 2p'{Q) \Vuj\^Q = 
[p{Q) + 2Qp'{Q)] |Vcu|' > K-\Q + ky |Vcu|' . (55) 
Thus (51) implies, via either (53) or (55) as appropriate, the inequality 

{uj, A {p{Q)uj)) < AH{Q) - K-\Q + k^ . (56) 

Applying eq. (46) to the left-hand side of (56) yields, for A = — {d5 + M) , 

(a;, A {p{Q)uj)) ^-*[uJ^*5d {p{Q)uj)] 
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= * [a; A d * d{f)ou)] = 

^_^y«(n+3) ^i^d[ujA *d{puj)] - [du A *d{pu)]} 

= [cu A *d{puj)] -*[vAljA *d{pu;)]} = n - Ts- (57) 

We express the first term in this difference, up to sign, as a divergence in the 
1-form dQ, writing 

Ti = *d[u) A *d{puj)] = 
*d[uj A* {p'{Q)dQ A u)] + *d[u A *pdw] = Tu + T12. (58) 
Equation (47) implies that 

T12 > ~|''"i2| = — \*d[u! A *{pv A cu)] I 

> -C {\Vuj\\v\p\uj\ + |V^;|pg+ |^;||cu||V(pc^)|) 

= C (-ri2i - ri22 - T123) ■ (59) 

We have, analogously to (54), the inequality p{Q) < K{Q + A;)*. Using this 
estimate and Young's inequality, we write 

- e|Vcup(Q + ky - c(£, k)\vWq + kyq. (eo) 

Kato's inequality and (48) yield, using \p'{Q) ■ Q\ < K{Q + ky, 

-T123 = -\v\\uj\\V{puj)\ = -\v\\uj\\p'{Q)VQ ■ uj + pVu;\ > 

-\v\\uj\ {\2p'{Q)\uj\V\uj\ ■ uj\ + \p{Q)Vu;\) > 
-2\v\\u;\ \p'{Q)-Q\ \V\u;\\ - \v\\u;\K{Q + ky \Vu;\ 
> -3\v\\u;\K{Q + ky\Vu;\ 
> -K{Q + ky (elVcup + Cie)\v\^Q) . (61) 
Substituting (60) and (61) into (59) yields, for a new s, 

T12 > -\ti2\ > 

- Ke{Q + ky\Vuj\^ - (C{s,K)\v\^ + K\Vv\) {Q + kyQ. (62) 
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Similarly, 

T2 — *[v /\UJ /\ *d{puj)\ > — C|i'||a;||V(pa;)|, 

which can be estimated by (61). Substituting (62) into (58), (58) into (57), 
and (57) into (56), and estimating T2 of (57) by (61) yields, again for a new 

*d[LuA* {p'{Q)dQ A iu)] ~ Ke{Q + A;)«| Vcjp 
-C(£, K){Q + ky {\Vv\ + \v\^) Q < AH{Q) - K-\Q + ky\S/uj\\ 
We obtain, choosing < £ < K''^ , 

< {K-^ -eK){Q + kf\Vu\^ 

< AH{Q) ± div (* [u; A * {p'{Q)dQ A uj)]) 

+C{Q + kf {^Vv\ + Q = L^(Q) + C(Q + k^ [\Vv\ + \v\^) Q. 

The ellipticity of the operator under condition (48) is obvious from 
the proof of Theorem 2. This completes the proof of Theorem 7. 

Corollary 8 Let {uj,v) be a solution of eqs. (46), (J^l) on where 
fl is a domain o/R", n > 2; E CC B CC fl is a compact singular set, com- 
pletely contained in a sufficiently small n-disc B, which is itself completely 
contained in Q. Let condition (48) hold. If n > A, let 2n/{n — 2) < p < n, 
where p is the codimension of S, and let u G L"'{B). If n = 3,4, let 
uj e L'^'^'>^{B)nL'^'^{B), where P = {p - e) / {p - 2 - e) for 2 < p < n, e > 0, 
and \<%<q. If {Q + kY{\Vv\ + l^l^) e L'^I'^{B) and iVvl + j-up e ^{B) 
for some p exceeding n/2, then uj is bounded on compact subdomains of Vl. 
This bound is uniform for k > 0. 

Proof. In (25) take $ = (Q + kY (|Vf | + \vf^ . Apply the arguments 
leading to (36) to show that Q is an H^'^ weak solution. Now choose^^ 

C = (K| + 5)^^-^^ 

for {uk} an increasing sequence chosen so that limk^^cuk = ui; f] & C^{B); 
r] > Q; 5 > Q; T > 1. Estimating (36) for this choice of test function implies 
in the limit that la;!"" e H^'^{B) for some r > 1. Also, (|a;|'^)'^ satisfies (25) 
for A < 2. Now Theorem 5.3.1 of Ref. 19 implies that is bounded. 
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7 The heat flow of solutions 

Consider the system 

-5{piQix,t))u;ix,t))^^^^, (63) 

du{x,t) = u{x,t), (64) 

where x G M, t G (0, T], and exterior differentiation is in the space directions 
only. Solutions of eqs. (63), (64) describe the heat flow, or gradient flow, of 
nonlinear Hodge maps. Notice that (64) implies du) = 0. 

If M is compact or if the normal component of u vanishes on dM, then 
the time decay of the energy 



is given by 



Et{uj) = - / p{s)dsdM 
2 Jm Jo 

' du \ , , , r I duj 



Equations (63), (64) imply that 

These identities together imply that 

^Et{u) = - I {dd{p{Q)u),p{Q)u) dM 
dt Jm 



{6 (piQp) , 6 {p{Q)uj)) dM < 0. (66) 

IM 

We conclude from (66) that a finite energy functional will remain so indefi- 
nitely. 

The local estimate for Q, taking M to be a bounded, open domain of R", 
is similar to its eUiptic counterparts in the proof of Theorem 7: If p {Q{x, t)) 
satisfies inequality (48), then 



{^p{Q) + Qp'iQ)) diQ^ - K-\Q + ky I V^l' > 
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(a;, A (pa;)) = *d[uj A* {p'{Q)dQ A a;)] + * a; A * 



du! 

'at 



using (63), (64), and (65), and 



IdQ 
2 dt 



di [[2P(Q) + Qp'iQ)) diQ\ ±dw*[ujA* {p'{Q)dQ A uj)] 




IdQ 
2 dt ' 



(67) 



This inequality is uniformly subparabolic whenever condition (48) is satisfied 
for /c > or g = 0. 

If M is a compact Riemannian manifold and -u : M x [0, T] — » iV, then we 
can obtain a global estimate for Q. In place of (14) we have the parabolic 
system 



Arguing as in the proof of Theorem 2, we add to the middle and right-hand 
side of eq. (20) a term of the form 



Let the sectional curvature of N be nonpositive. We obtain as in (67) the 
inequality 



where Cr depends on the Ricci curvature of M. This inequality is, of course, 
also uniformly subparabolic whenever condition (17) is satisfied. In fact 
we have an a priori estimate in this case, which strongly depends on the 
ellipticity of L^^. 

Theorem 9 Let u{x,t) be a mapping of a smooth, compact n- dimensional 
Riemannian cylinder M x [0, T] into N, where N is a smooth, compact Rie- 
mannian manifold of nonpositive sectional curvature and T is a finite num- 
ber. Suppose that u smoothly satisfies (68) with p small in L'*/^(M), s > n, 
and with p'{s) < 0, s G [0,(5]- Let condition (17) hold for each point 
(x, t) e M X [0, T] and let Q{x, 0) < 1. Let the H^'^ Sobolev inequality hold on 





luQ) + CrpQ---^>o, 



(68) 



28 



M for constants Si, 82- Then there is a constant c{s, K, q, M, N, T, Si, S2) 
such that for q > 0, 



sup (sup Q{x, t)] < cr"/2(''+^) {E [uj{x, 0)]) 

te(0,T] \x€M J 



1/(9+1) 



where E is the nonlinear Hodge energy. 



Proof. Multiply inequality (68) by {Q + ^ for r > 1 and /3 > 0. 
Replace the time derivative in (68) by the (identical) time derivative of + 
and integrate over M. We obtain 



^"'^ / iQ + PydM<[ {Q + I3y-^V ■{a{u:)VQ) dM 
ut Jm Jm 

+ C f p-{Q + /3Y dM, 

JM 



(69) 



where V is the gradient on M and a is the matrix- valued function of inequal- 
ity (27). Because M is compact, Stokes' Theorem implies that 

/ (g + py-^ V ■ (a {uj) VQ) dM = f V-(a{uj){Q + PY'^ VQ) dM 
Jm Jm ^ ' 

-[ v({Q + PY-^)-a{u)VQdM^- [ {r - 1) {Q + PY'^ a{u)\VQ\'' dM 
Jm ^ ' Jm 

<-mi[ {r-l){Q + PY~^\^Qf dM^-mi [ V ({Q + (5Y~^)-VQ dM, 
Jm Jm ^ ' 

where mi depends on k, K, and q. Now 



v((Q + /3) 



r/2 



^(g + /3)(^-^)/^vQ 



-^{Q + (iY-\^QY 



4 (r - 1) 

so we can write inequality (69) in the form 



— -V ((g + /5)^-^) vg, 



dt 



f {Q + PY 

Jm 



dM < 



4mi(r — 1) 



/Jv(«5 + «"^) 



dM 



+C\\pl/2\\{Q + PY\\l/is-2)- 
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Employing the parabolic DeGiorgi-Nash-Moser iteration as in Sec. 4 of Ref. 
25, taking pq — q + 1, we obtain, letting (3 tend to zero. 



sup (supg(a;,t) ) < Cr"/2(9+i) /■ \Q{x,0)\''+Um 



te{o,T] \xeM 
Because 



^ 1/(9+1) 



^ {sp{s)) = p{s) + sp'{s) < p{s), 



we have, for p'{s) < 0, the inequality 



Qp{Q) = l ^M^))^^^/ ^(^)^^- 



Thus 



'^E\t=Q> / / ' p(s)dsdM> / Q(x,0)p(Q,0)dM > 
Jm Jo jm 

[ Q{x, 0) [p{Q, 0) + 2Q{x, Q)p'{Q)] dM > R-^ [ Q{x, oy+^dM. 
Jm Jm 

Taking the (g + 1)** root of this inequahty and using (66) completes the proof 
of Theorem 9. 

A local version of Theorem 9 would argue from inequahty (67) rather than 
(68). The initial argument is as in the proof of Theorem 9 except that the 
integration is against cut-off functions. The Moser iteration is implemented 
as in Ref. 26. 
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